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Abstract 

A recently proposed method for a strong coupling analysis of scattering phenomena in 
hamiltonian lattice field theories is applied to the SU(2) Yang-Mills model in (2 -|- 1) 
dimensions. The calculation is performed up to second order in the hopping parameter. 
All relevant quantities that characterize the collision between the lightest glueballs in the 
elastic region - cross section, phase shifts, resonance parameters - are determined. 



Introduction 



This paper is a continuation of ref. [1] where I have developed a method to obtain strong 
couphng expansions for scattering quantities in hamiltonian lattice field theories. The keystone 
of the formalism is the derivation of an effective quantum mechanical system using Bloch's 
perturbation theory [2]. This system determines the two-particle dynamics of the theory by 
means of a free hamiltonian and an interaction potential. Thereby the problem of scattering in 
lattice field theories is reduced to a quantum mechanical problem of scattering at a potential 
and can be solved using the powerful methods of ordinary quantum mechanics. 

In order to improve our understanding of the low-energy regime of elementary particle 
physics several authors have studied collision phenomena in lattice field theories in recent years 
[3] - [12]. It is the intent of my approach to supplement these investigations by a discussion 
within the framework of strong coupling perturbation theory. 

In [I] the [d-\- 1) dimensional Ising model was chosen as a test case. The objective here is to 
investigate the collision of glueballs in non-abelian gauge theories. These models are described 
by the famous Kogut-Susskind hamiltonian [13], and a strong coupling perturbative analysis 
was first carried out by Kogut, Sinclair and Susskind in their pioneering work on the SU(3) 
glueball spectrum [14]. It is characteristic for Yang-Mills theories that the spectrum is rather 
extensive (to get an idea of the variety of Wilson loops the reader is referred to ref. [15]). 
In the static limit the excitations are highly degenerate. In particular, one finds a number of 
single-particle states that have the same static energy as two-glueball configurations. When the 
next-to-nearest neighbour interaction is turned on, this sector is deformed into a two-channel 
system describing two identical particles coupled to an additional channel of one-particle states. 
The reduced hamiltonian is a 2 X 2 matrix in the space of the two channels. 

As far as the strong coupling analysis of glueball scattering is concerned this observation 
adds a new aspect compared to the Ising model. In the present paper I will demonstrate to 
what extent the methods developed with the Ising model must be broadened to be applied 
to an investigation of non-abelian gauge theories. It will be argued that the effect of the 
coupling to the single-particle channel can be taken into account by an additional potential in 
the two-glueball channel. 

Since there are a lot more different Wilson loops in three than in two spatial dimensions 
I decided to discuss the conceptual ideas with the SU(2) gauge theory in (2 -|- I) dimensions. 
The present paper is not self-contained in the sense that the reader is referred to [I] from time 
to time. 

In section I the model is introduced. I discuss the static limit and analyse the dynamics 
of the lightest glueball up to second order in the hopping parameter k. Scattering between 
the lightest glueballs is supposed to take place in the sector of eight-link excitations. Here one 
distinguishes between the connected Wilson loops and the disconnected configurations. The 
former are interpreted as single-glueball states whereas the latter represent two distant light 
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particles. 

Section 2 is the heart of this work. By an iteration of Bloch's method I show that the eight- 
hnk sector decomposes into a number of bound states and an effective two-channel system where 
the light glueballs are coupled to three heavy particles. As far as the two-glueball channel is 
concerned an effective Lippmann-Schwinger equation is derived and a series representation for 
the scattering solutions can be found. The leading order approximation for the scattering 
amplitude is calculated explicitly. 

In section 3 the analysis is completed by the computation of the phase shifts for the dif- 
ferent symmetry sectors. The techniques developed with the Ising model require only a slight 
modification. A resonance is detected in the collision of the light glueballs and discussed in 
some detail. 

I end with some concluding remarks and an outlook. In order to make the paper more 
readable I have collected some supplementary material as well as a few merely technical items 
in three appendices. 
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1 SU(2) gauge theory in (2 + 1) dimensions 



1.1 Preliminaries 

We consider the SU(2) Yang-Mills theory on a rectangular lattice 

T = {x\x e7L\ -L/2 <xk< L/2, k = 1,2} (1.1) 

with linear extent L and lattice spacing a = 1. The directed links between neighbouring sites 
are labeled by a site x and a spatial index k such that the sites x -\- k and x locate the 
beginning and the end of the link (a;, k). By definition, (a;, k) and [x + A;, —k) denote different 
links although they occupy the same lattice bond. A classical SU(2) lattice gauge field is an 
assignment of a link variable U{x^ k) G SU(2) to any bond (a;, k) of the lattice V such that 

U{x,k)-^ = U{x + k,-k) . (1.2) 

On a finite lattice one has to specify boundary conditions. For definiteness we choose them to 
be periodic, hence 

U{x,k) = U{x + Li,k) , /=1,2. (1.3) 

In the quantum theory the classical field U{x^ k) becomes an operator field U(a;, k) acting 
on some Hilbert space Ti' of states. At each link (a;, k) of the lattice V we associate a copy of 
iv^(SU(2)), the space of square integrable functions on the group manifold (cf. appendix A), and 
call it iv^(a;, k). For a finite lattice the Hilbert space Ti' is defined as the tensor product of all 
L^[x^k). In this Schrodinger picture the states are complex wave functions f[U] that depend 
on the basic field variables, i.e. the argument runs over all (periodic) gauge fields U{x^ k) on F. 
The group of (periodic) gauge transformations g{x) on F acts on the wave functions through 

f[u] fm , (1.4) 

where U^{x, k) denotes the gauge transform of the link field 

U'{x, k) = g{x) U{x, k) g{x + k)-^ . (1.5) 

The canonical field operators are identified as multiplication operators on 7i' 

[l]{x,kUf][U] = U{x,kUf[U] . (1.6) 

In appendix A we introduce the Casimir operator F^ of the gauge group SU(2). By F^(a;, A;) 
we denote a copy of F^ which acts on L'^{x, k) as F^ and as the identity operator on the other 
components of the tensor product. Now we are prepared to write down the pure SU(2) lattice 
gauge theory hamiltonian which was first derived by Kogut and Susskind [13] 
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Here we have abbreviated the plaquette operator 

P(a;) = Tr [u(a;,l)U(a; + i,2)U(a; + 2,l)-^U(a;,2)-^] . (1.8) 

and the so-called hopping parameter n = In addition there is the discretized version ol 

Gauss' law with zero external charge. By this constraint the physical Hilbert space is identified 
with the subspace of 7i' that only contains the gauge invariant wave functions. We denote the 
physical Hilbert space by 7i. 

Let us have a closer look at the structure of the Kogut-Susskind hamiltonian. The first part 
of (1.7) is a sum over local operators that act non-trivially on one component of the tensor 
product 7i' only 

IHo = EEE2(«^'^) • (1-9) 

X k=l 

Therefore it is a static term whereas 

\u, = -j2nx) (1.10) 

X 

couples the links of the plaquettes and provides a kinetic operator. So we are back with the 
situation of the Ising model and ready for a perturbative investigation in the hopping parameter 
K. For convenience we introduce 

IH = (2//)IH' = IHo + kIHi . (1.11) 

1.2 The static limit 

The static hamiltonian is a sum of single-link Casimir operators. From the analysis in appendix 
A we conclude that the eigenfunctions of IHq are products of representation functions of the 
gauge group SU(2). Clearly speaking, for each link (x^k) on the lattice F one chooses a 
quantum number n that denotes the dimension of the irreducible representation and indices 
a^b = 1, . . . n. Then one forms the product of the associated representation functions [/"^(a;, k). 
The eigenvalue of this static eigenfunction is the sum of the corresponding one-link eigenvalues 
(l/4)(n2-l). 

The ground state of IHq is obtained if we set n = 1 for all links, i.e. we choose the trivial 
representation at each link. The corresponding wave function f[U] = 1 is gauge invariant, so 
we are led to interpret this state as the static, physical vacuum and denote it by 

IHo|O)° = 0, °(0|0)° = 1. (1.12) 

The low-lying eigenstates of IHq have few excited links. Let the operator corresponding to the 
representation function [/"^(a;, k) be denoted by a Roman letter U"5(a;, k). Then the static state 
with one excited link (a;, k) is defined as 

\x,k;n,ab)'' = N V^bix, k)\n)'' . (1.13) 
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Figure 1.1: Static one-particle excitations. 



To determine the value of the normalizing constant N we have to calculate 

°(a;, k;n,ab\x, k; n, a6)° = iV' °(0 | U^Jx, hf \^l^,{x, k) \ 0)° . (1.14) 

The method to evaluate such matrix elements relies on the Clebsch-Gordon series for the 
product of representation functions and is described in detail in appendix A. Here we obtain 
N = \/n. From the definition (1.13) it follows 

IHo|a;, A;;n,a6)° = (l/4)(n2 - \)\x,k;n,ahf , (1.15) 

i.e. for each link with quantum number n > 1 the colour electric energy is increased by the 
amount (l/4)(n^ ~ !)• We say that a colour electric fiux has been created, and the quantum 
number n is a measure for the strength of this fiux. 

Within this picture the requirement of gauge invariance for the physical states simply means 
that the colour electric fiux out of any site of the lattice F is zero. Therefore we have to form 
appropriate products of the non gauge invariant states (1.13) to generate the physical spectrum 
of the static hamiltonian IHq. For instance, a gauge invariant excitation is built by taking the 
trace of n = 2 fiux elements around one single plaquette 

|a;)° = V{x)\^f , (1.16) 

where the plaquette operator was defined in eq. (1.8). These states have energy 4°) = 4- (3/4) = 
3. They satisfy 

""{x' \xf = b{x' -x) . (1.17) 

and are interpreted as the lightest static glueballs. Since the trace of SU(2) matrices is real 
the fiux along the plaquette does not have an orientation. The different one-particle states are 
distinguished by the vector x indicating the position of the plaquette on the lattice. 

Heavier glueballs of mass 4°) = 6- (3/4) = 9/2 are created from the vacuum by six-link loop 
operators in the same way as the four-link states (1.16). It is most suggestive and convenient to 
define them pictorially as in fig. 1.1. The two states |a;, 1)° and |a;, 2)° have the same geometrical 
shape and are transformed into each other by lattice rotations G 0(2, TL). 

The next glueballs in the static spectrum are the loops of length eight with energy E^^^ = 6. 
There exist five geometrically distinct classes which are invariant under lattice rotations. We 
adopt the notation 

\x,{k,ak))\ (1.18) 

where A; = 1, ... ,5 indicates the class and ak distinguishes between the different states of the 
same shape. In fig. 1.1 we depict the five prototypes. The number in brackets denotes the 
dimensionality of the class. 

Since the colour electric energy of the eight-link glueballs is just twice the energy of the 
four-link states the static eigenspace of energy 6 also contains two-particle excitations 

\x\xy = V{x^)V{x^)\n)'' , \x^-x^\>l. (1.19) 
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The above restriction is necessary to avoid that the two excited plaquettes overlap on one 
or more hnks. In case they do products of representation functions occur which have to be 
reduced by means of the Clebsch-Gordon series. The resulting states will be closed loops which 
carry different flux on the (former) common links and represent single-glueball configurations 
belonging to other static eigenspaces. Since the plaquette operators at different sites commute 
it follows 

\x\xy = \x\xy (1.20) 

and the normalization is 

'{y\y' I x\xy = '{y^ \ x^'iy' \ x')' + '{y' I ^Tiv' I ■ (1-21) 

Above the two-particle threshold there are seven-link excitations |a;, k; 7)° whose prototype 
is depicted in fig. f.f. Such configurations are created by acting with neighbouring plaquette 
operators ^^{x) and V{x -\- k) on the vacuum and taking the n = 3 representation on the 
common link, as indicated by the dashed line. The static colour electric energy is Ef"^ = 
6 • (3/4) + 2 = 13/2. 

This overview illustrates the richness of the (static) spectrum of lattice Yang-Mills theories 
already in two spatial dimensions. If we go to higher energies we will encounter excited glueball 
states (n > 2), multi-particle configurations and rather long n = 2 fiux loops. 

It is the aim of the present work to investigate the elastic scattering between the lightest 
glueballs (f.f6) by means of a perturbative expansion up to second order in the hopping pa- 
rameter K. As a prerequisite one has to analyse the one-particle dynamics of these excitations. 
Thereby I repeat the strong coupling technique for gauge theories as described in ref. [14]. 

1.3 Vacuum sector 

Since the physical quantities are only defined relative to the ground state of the theory we first 
have to calculate order by order in k, the shift of the static vacuum when the colour magnetic 
perturbation (f .10) is turned on. From the general formula (2.14) in [I] we obtain the reduced 
hamiltonian for the vacuum sector 

= /cPolH^Po + ^'PolHi^^rf^lHiPo + 0{k') , (1.22) 

— IHo 

where the projection operator is given by Pq = |0)°°(0|. The action of the plaquette operator 
^(a;) on the vacuum creates an elementary four-link configuration which is orthogonal to the 
vacuum eigenspace. Consequently there is no first order contribution to Hq. 

As for the Ising model we introduce a pictorial description of the processes in strong coupling 
perturbation theory. In fig. 1.2 the prototype of the non- vanishing matrix elements contributing 
to the reduced hamiltonian (1.22) in second order is visualized by going from the left to the 
right. The initial state is the static vacuum without any colour electric fiux. The action of 
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Figure 1.2: Second order vacuum process. 

the plaquette operator 'P{x) is represented by a right-arrow and leads to an intermediate four- 
link configuration of energy 3 at site x. Applying 'P{x) once more the intermediate state is 
annihilated and we end up with the vacuum as the final state. The value of each such matrix 
element combined with the energy denominator is —1/3. Because there are squares on a 
two-dimensional periodic lattice of linear extent L the counting factor in . One obtains 

H^ = -«2(x2/3)Po + (9(/c3) (1.23) 

as the effective hamiltonian and 

E, = -K^ {L'/3) + 0{^) 

for the ground state energy. 

1.4 The lightest glueball 

The projection operator on the static four-link space S4 reads 

P4 = EI«^)°°(«^I • (1-25) 

X 

From the general Clebsch-Gordon series (A. 21) we know that 2(x)2 = 1®3, i.e. the product of two 
fundamental representations does not contain the fundamental representation itself. Therefore 
the action of the plaquette operator on a four-link state leads out of the static subspace S4. As 
for the vacuum sector there is no first order contribution to the reduced hamiltonian. 

To second order in perturbation theory one encounters five different types of transitions. 
The corresponding graphs are collected in fig. 1.3. Column (a) represents the process where the 
initial state remains undisturbed while a vacuum fiuctuation occurs on a distant square. The 
value of each such matrix element is equal to the vacuum process depicted in fig. 1.2, but the 
counting factor is reduced by the number of plaquettes that have one link in common with the 
initial state. In fig. 1.3 (b) the complementary graph is shown where the vacuum fiuctuation 
now has one link in common with the initial state. This link may either be in the trivial 
representation or in the adjoint representation leading to an intermediate six-link configuration 
|a;, or a seven-link excitation |a;, k; 7)° (cf. fig. 1.1). The author hopes that the diagrammatic 
description of the different processes is for the most part self-explanatory such that we must 
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Figure 1.3: Second order contributions to the static four-link sector. 
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not go through all of the graphs in fig. 1.3. In appendix A all types of elementary matrix 
elements which are needed to first and second order in perturbation theory throughout this 
thesis are given explicitly. With these results at hand the evaluation of the different diagrams 
is straightforward. To be as transparent as possible the values of the matrix elements combined 
with the energy denominator as well as the counting factors for the processes are also listed in 
fig. 1.3. 

We deduce the reduced hamiltonian 

h; = [-(ivV3) + 29/105] P4 - (KV21) E E (l« + W + 1* - W) + 0{k^) ■ (1-26) 

X k=l 

The first summand is the result of the graphs fig. 1.3 (a) - (c) where the final state is equal 
to the initial one. The processes fig. 1.3 (d) - (e) result in a shift of the initial configuration 
by one plaquette in one of the four spatial directions and lead to the second term of (1.26). 
Subtracting the vacuum energy shift we define the infinite volume hamiltonian H4 = H4 — Eq 
and pass to the coordinate space representation 

H4 = (3/35) + (KV2I) A + 0{k^) , (1.27) 

where the lattice laplacian was defined in [1, eq. 2.29]. The eigenfunctions are plane waves 
exp[ipx) and the dispersion relation for the lightest SU(2) glueball in two space dimensions 
reads 

E^{p) = 3 + K^ (3/35 + pV21) + 0{k^) . (1.28) 
For the physical rest mass and the kinetic mass of the particle we obtain 

m^ = 3 + {3/35)K^ + 0{k^) , m^^ = [2/21) + 0{k^) . (1.29) 

The reduced hamiltonians for the six- and seven-link sector are obtained in much the same 
way. In particular, it is straightforward to see that there are no first order contributions either. 
To second order one finds that the non- vanishing matrix elements (combined with the energy 
denominator) corresponding to the processes that either shift the initial configuration by one 
plaquette or transform between the two different states in Sq and 87^ respectively, add up to 
zero. Consequently, the heavier particles remain static and degenerate. The mass is 

mr,6 = 9/2 + (8/21) + 0(^3^ ^^_3Q^ 

for the six-link states and 

mr,7 = 13/2 + (8/21) + 0(^3^ ^^_3^^ 
for the seven-link configurations. 
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2 Scattering of glueballs 



2.1 The eight-link sector as a two-channel system 

The static subspace Eg, of energy 6 contains the connected loops of length eight depicted in 
fig. 1.1 together with pairs of separated four-link loops defined through eq. (1.19). As for the 
Ising model I restrict the analysis to the center of mass system. Then the two-particle states 
depend on the relative coordinate r = — with |r| > 1, and the properties (1.20) and 
(1.21) translate to 

|r)° = |-r)°, %r' \ = 8{r' - r) + 8{r' + r) . (2.1) 

In view of the forthcoming calculations it is advantageous to distinguish between the configu- 
rations with a relative distance |r| > 2 and those depicted in fig. 2.1. The connected eight-link 
states are independent of any lattice vector in the center of mass system and are denoted by 
|A;,a/;)°. To fix the notation they are listed in fig. 2.2. 

Before we can write down a basis for the linear space Eg, we have to remember a specialty in 
SU(N) gauge theories, the Mandelstam constraints [16]. These are non-linear identities among 
loops which intersect or touch at at least one point. Their origin are basic properties of SU(N) 
matrices. In case of = 2, all Mandelstam constraints can be derived from the following 
fundamental identity for arbitrary SU(2) matrices U and V 

Tr[/y + Tr[/yt-Tr[/Try = . (2.2) 

It is an easy exercise to verify that (2.2) leads to linear relations between certain states in Eg, 
namely 

|4,l)° + |5,l)°-|i + 2)° = , |4,2)° + |5,2)°-|i-2)° = . (2.3) 

We agree to choose {|A;,a/;)° | A; = 1 . . . 4} together with the two-loop configurations fig. 2.1 as 
linear independent states. However, they are not orthogonal, rather 

°(4, 1 I i + 2)° = °(4, 2 I i - 2)° = 1/2 . (2.4) 

In order to obtain an orthonormal basis for the eight-link space we define 

|4', 1)° = (V3/3) (2 |4, 1)° - |i + 2)°) , |4', 2)° = (V3/3) (2 |4, 2)° - |i - 2)°) . (2.5) 

For convenience I introduce a new notation for the states fig. 2.1 

\h',kf = \2-k)\ |6',1)° = |i + 2)° , |6',2)° = |i -2)V (2.6) 

Below the prime will be dropped whenever there is no danger of confusion. 
We define the operators 

p:c = (i/2) E i^r(^h K = i:i:\^.a,r{Ka,\. (2.7) 

|r|>2 k=l a.k 
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|2-1)° |2-2)° |l + 2)° |l-2)f 



Figure 2.1: Two-glueball configurations with a small relative distance 1 < |r| < 2. 

The first one, P^^, projects onto the static two-loop configurations with a relative distance larger 
than 2. When the kinetic hamiltonian kIHi is turned on, they are deformed into states describing 
the relative motion of the light glueballs. On the contrary, the states a/;)° | k = 1, . . . , 6} 
stay at rest in the center of mass system and are interpreted as single-particle excitations. 
Thereby one is led to a natural decomposition of the Hilbert space 

Ss = ® K , (2.8) 

where Ti'^^ and denote the orthogonal subspaces associated to the above projection operators. 
Accordingly, the states of the system are represented as two-component vectors 

Ix) = j , IV-) e K. , W) e K ■ (2.9) 

The infinite volume reduced hamiltonian is a 2 X 2 matrix 

^ H' V 

sc cc 

V' ^ M' 



(2.10) 



v 



cc 



Wl 



ith 



C = <clHiP:c + «'P:clHi^^lHiPL-i?oPL + 0(«^), (2.11) 

M' = «p;,IH,PU«XlHi^^lHiP;,-i?oPU0(«') , (2.12) 

= «P:JH,PU«'p:cIHi^^IHiPU0(«') . (2.13) 
Here Ps represents the projection operator on the entire Hilbert space Eg,, hence 

P8 = p:c+Pb. (2.14) 

12 



1)° 



|2,1)° 



|2,2)° 



|2,3)° 



|2,4)° 



|3,1)° 



|3,2)° 



|4,1)° 





|4,2)° 



|5,1)° 



|5,2)° 



Figure 2.2: Single-glueball configurations of static energy 6. 
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To summarize, we are concerned with a two-channel system that describes two identical particles 
coupled to an additional channel containing 13 single excitations. The hamiltonian H^^ is the 
operator of the kinetic energy of the two particles possibly supplemented by an interaction 
potential. In the single-particle channel the hamiltonian M' is just a mass operator. Finally 
the operator V^^ provides the coupling between the two channels. 

2.2 The reduced hamiltonian 

This section is devoted to the deduction of the reduced hamiltonian (2.10). We start with the 
two-particle channel. 

The first order contributions to H^^ come from processes where the action of the plaquette 
operator ^(a;) on a two-glueball initial state leads to a two-glueball final state. From our 
considerations in section 1 we know that such a transition is not possible. 

To second order the perturbation acts on either of the two glueballs as described by the 
graphs in fig. 1.3, whereas the other four-link loop remains undisturbed. The only thing we 
have to take care of is the Bose symmetry of the two-particle configurations. This is guaranteed 
by a four times larger counting factor for the graphs corresponding to fig. 1.3 (b) - (e). One 
obtains 




+ ^^' E [(|U/)° + |2-^)°)°(2 4 + /| + (|^-/)° + |2 J)°)°(2 4-/|] 

+OiK^) . (2.15) 

Since is defined on the subspace of two-glueball configurations with a relative distance larger 
than 2 we have to subtract those contributions to the sum over |r| > 2 that lead out of 
This is the origin of the additional terms in (2.15). 

Next we turn to the mass operator (2.12) in the single-particle channel. Unlike for the two- 
particle channel first order transitions exist and are shown in fig. 2.3. The plaquette operator 
creates a four-link excitation on the lower left square having two links in common with the 
initial state |2, 1)°. Since the product of two fundamental representations contains the trivial 
one we obtain another eight-link configuration, |4,2)° or |6,2)°, as the final state. Using the 
methods of appendix A the value of the corresponding matrix elements are computed to be 1/2 
for the graph fig. 2.3 (a) and 1/4 for the graph fig. 2.3 (b), respectively. Further transitions 
between states in 7i[^ are |1, 1)° — )• |2, 02)° belonging to the class fig. 2.3 (a) and |3, — )• |5, 
whose prototype is fig. 2.3 (b). Adding up all contributions the first order reduced hamiltonian 
in the single-particle channel reads 

M; = K [M^ + M^^) , (2.16) 
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Figure 2.3: First order transitions. 



where we have abbreviated 



^1 = -If2x:i2,^r(i,ii+V3x:i2,^r(4,2i+V3x:i2,^r(4,ii 
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The second order contributions to (2.12) can be grouped together in seven classes, repre- 
sented by the graphs fig. 1.3 (a) - (c) and fig. 2.4. The values lor the corresponding matrix 
elements are the same lor all different processes belonging to the same class whereas the count- 
ing lactors also depend on the shape ol the initial and final states. We get 



M'2 = (M^ + M^^) + 0{k 



(2.17) 



with 
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Figure 2.4: Second order corrections to the static eight-link sector. 
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Finally we have to calculate the operator (2.13) that determines the coupling between the 
two channels. To second order in n the single-particle excitations {|A;,a/;)° | A; = 4,5,6} are 
coupled to the two-glueball configurations {\rY | 2 < |r| < 3} via the processes fig. 1.3 (d) - 
(e) that shift one of the excited plaquettes by one square. The result reads 

VL = -|r-^E|3-^n5,^|-|-.^ E (|2.U/)° + |2.^-/)°)°(5,^| 
k=i '^^ k^i=i 

+ 1^^'^' E [|2-U/)°(°(4,1|-V3°(6,1|) 
""^ k^i=i 

+ |2 4-/)°(°(4,2|- V3°(6,2|)] +0{^) . (2.18) 
2.3 Diagonalization of the reduced hamiltonian 

In the preceding section we have deduced the second order approximation of an operator that 
determines the dynamics in the eight-link sector Sg- Perturbatively, the eigenvalues and eigen- 
states of this reduced hamiltonian can be obtained by an iteration of Bloch's method. Clearly 
speaking, one diagonalizes the first order part exactly and regards the higher order contribu- 
tions as a perturbation. In case that there are eigenspaces of different energy the problem is 
further reduced leading to new effective hamiltonians that act on lower dimensional Hilbert 
spaces. 

The first order contribution to Hg is confined to the single-particle channel and represented 
by the mass operator (2.16). Among the 13 states {|A;,a/;)° | k = 1,...,6} the degeneracy 
is partly lifted as 7i[^ decomposes in seven orthogonal subspaces. To begin with one finds 
four one-dimensional and two two-dimensional eigenspaces with masses different from 6. So 
these excitations decouple from the two-particle channel and represent bound states of the 
system. Regarding M2 as a perturbation the second order correction to the various subspaces is 
determined straightforwardly by means of corresponding reduced hamiltonians. The degeneracy 
is not reduced any further, and the associated masses are 





= 6- 




= 6 - 1.225 /€ + 0.676 + 0(/€3) , 


(2.19) 


me 


= 6- 


f . + |i.^ + (9(.^) 


= 6 - 0.707 K + 0.487 + OiK^^) , 


(2.20) 


nir 


= 6- 


1 . + ii + (9(.^) 


= 6 - 0.25 K + 0.523 + 0{k^) , 


(2.21) 


nig 


= 6 + 


1 . + ii + (9(.^) 


= 6 + 0.25 /€ + 0.523 + 0(/€^) , 


(2.22) 


nig 


= 6 + 


f . + li + (9(.^) 


= 6 + 0.707 K + 0.487 + OiK^^) , 


(2.23) 


mio 


= 6 + 


f . + ^.^ + (9(.^) 


= 6 + 1.225 K + 0.676 + 0{k^) . 


(2.24) 



So for this part of the spectrum the job is done. 
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In addition, there is a five- dimensional space Ti^,' of single-particle excitations that remain 
degenerate with the two-glueball channel to first order. One obtains a modified two-channel 
system where the number of additional states coupled to Ti'^^ is reduced from 13 to 5. In order 
to construct a modified reduced hamiltonian that acts on the reduced space of the two channels 
® '^h ^iid thus determines the scattering of the two four-link glueballs we specify a basis 
of normalized single-particle states 



= 

|2') = 
|3') = 
|4') = 
|5') = 

Introducing the projection operator 



1/2)(|4,1)°-V3|6,1)°) , 
1/2)(|4,2)°-V3|6,2)°) , 

2 2 

V3/3) |1, 1)° - (1/2) |4, - (V3/6) |6, kf 

k=i 

V2/2) (|2,1)°-|2,2)°) , 
V2/2) (|2,3)°-|2,4)°) . 



k=i 



K = j:\k'){k'\ 



(2.25) 
(2.26) 
(2.27) 
(2.28) 
(2.29) 

(2.30) 



k=i 



the desired reduced hamiltonian is given by 

/ 



p/ jj/ p/ 

sc sc sc 

p// y/ t p/ 

b cc sc 



p/ y/ p// 

SC CC b 



(2.311 



From (2.15) it follows that H^^ = Pg^Hg^P^^, so nothing changes for the two-particle channel. 
By means of the explicit formula (2.17) the mass operator for the reduced single-particle channel 
is easily determined 

p"M',p;: = iSi'^^ Ei^')(^'i + ii'^^|3')(3'i-ii-^E(i^')(3'i + |3')(^'i) 



k=i 

5 



k=i 



+m-'j:w){k'\+o{K' 

k=4 



(2.32) 



Having a closer look at the channel coupling potential (2.18) one finds that the first order states 
(2.25) and (2.26) appear explicitly. Hence, for the coupling between 7i'^^ and Ti'-^ one gets 

KK.K = ^V3'^' E (|2J + /)°(l'| + |2.^-/)°(2'|) . (2.33) 
""^ k^i=i 

As far as the single-glueball channel is concerned the diagonalization of the operator (2.32) 
leads to four subspaces of different mass 

6 + (fi - "^HF) + 0{k') = 6 + 0.449 + 0{^) , (2.34) 
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nil 



= 6+(f|i + ^^i||p)/c2 + (9(/c3) = 6 + 1.216/c2 + (9(/c3), (2.35) 

ms = 6 + f|i| + 0(/€3) = 6 + 1.175K2 + 0(/€3) , (2.36) 

m4 = 6 + i|| + 0(/€^) = 6 + 0.673K2 + 0(/€^) . (2.37) 
The first three are of dimension one and can be represented by the normaiized states 

|1) = /3i (|l') + |2') + «i|3')) +C(/€2) , (2.38) 

|2) = /32 (|1') + |2') + «2|3')) +C(/€2) , (2.39) 

|3) = - (V2/2) (if) - |2')) + 0{k') . (2.40) 
The aigebraic coefficients have the vaiues 



^, _ 3160 I V12552978 _ c; Qi R ^, _ 3160 V12552978 _ n QQS 

"1 - 1133 + 1133 - ^-^^^ ' "2 - n33 113^ " -^-^JO 

and 

9 \-l/2 ^ / 9 \-l/2 



/3i = (^a^ + 2j = 0.164 , (3^ = [ai + 2) = 0.688 . 

The remaining subspace, associated with the mass m4, is of dimension two and spanned by the 
vectors 

|4) = \i') + Oi^) , |5) = \5') + Oi^) . (2.41) 

From the expiicit formuia (2.33) it foiiows that the states (2.41) do not coupie to the two-giuebaii 
channei They represent two further bound states of the system. 

Thus, to second order in the hopping parameter k, one is left with an effective quantum 
mechanical system which describes two identical light glueballs of mass nij- = 3+0.085 k,'^ + 0{k,^) 
coupled to an additional channel that contains three single, heavy particles (2.38) - (2.40) of 
mass nii > 2my. We define the projection operator 

Php = El^)(^l (2.42) 
k=i 

and denote the heavy particle channel by ?ihp- For notational purposes it is convenient to 
express the masses of the heavy particles as follows 

m, = 6 + /€V.- for z = 1,2,3. (2.43) 

Then the effective hamiltonian that determines the dynamics in the space of the two channels 
^sc ® ^hp reads 

Hs= I , (2.44) 



VJ M 
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where the mass operator in the heavy particle channel is given by 

k'M = Pi^p V'i n Pi^p = K'j2l^k \k) {k\ , (2.45) 

k=i 

and the coupling between the two channels is mediated through the potential 

V,, = P;: Php = VL + 0{k^) (2.46) 

with the leading order contribution 

VL = E (|2-U/)° + |2 4-/)°) (/3i(l|+/32(2|) 

-Ij^ E (|2 4 + /)° - |2 4 - 0°) (3| . (2.47) 

Note that the coefficients in eq. (2.47) are about 0.08, i.e. the strength of the coupling between 
the two channels is comparatively weak. 

Compared to the previous result (2.10) the effective operator (2.44) is obviously less com- 
plicated. So the reduction of the dimensionality really simplifies the problem. 



2.4 Lippmann-Schwinger equation 

As a prerequisite for a determination of the scattering solutions in the coupled two-channel 
system described above one has to identify the free dynamics. Concerning the two-particle 
channel we follow the strategy already applied to the Ising model. Any state G 'H'^^ has a 
unique expansion 

IV') = ^EV'(^)k)°, (2.48) 

where the complex lattice function V'(^) is symmetric and subject to the restrictions V'(^) = 
for |r| < 2. The space of these (normalizable) functions is identified as the coordinate 
representation of the Hilbert space 'H'^^. As emphasized in section 2.5 of [1], a wave function 
that vanishes at certain lattice points near the origin cannot describe the independent relative 
motion of two particles. Instead we consider the space of unrestricted wave functions Tisc — 
L'^{T1?\ For technical reasons we have also dropped the restriction for the wave functions to be 
symmetric. At the end of our calculation the true two-glueball sector is restored by projecting 
on the subspace of symmetric states (cf. the detailed discussion with the Ising model in [1]). 

In order to construct a self-adjoint extension of the reduced hamiltonian (2.15) on Tigc we 
introduce the operator 

Dnc = Do + E (d^ + + D,.^ + J + + D^.^ + D.^^^ + D.^.^ . (2.49) 
fc=i 
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It projects on the wave functions that vanish outside a disc of radius 2 around the origin. The 
extended hamiltonian is defined as 

H,, = (]l - C (]l - + D(A) , (2.50) 

with the local multiplication operator 

D(A) = AoDo + X: [Afc (d^ + D_^) + A,+2 (d^.^ + D_2.^ 

k=i 



+ A5 + 2 + D-i-2j + ^6 (Di_2 + D_i + 2 

and arbitrary real constants A^, // = 0, . . . , 6. The explicit result is deduced straightforwardly. 
One finds 

H,c = /c' (T + V,c) , (2.51) 

with a kinetic term 

T = 6/35 + (2/21) A + 0(/€) , (2.52) 

and a short range potential 

V.C = -Dnc T - T D„c + D„c T D„c + D(A) + 0{k) . (2.53) 

As the reader will have expected the kinetic part is just twice the effective hamiltonian (1.27) 
for the four-link sector. Therefore (2.52) describes the free propagation of the light particles 
in the two-particle channel. The corresponding (improper) eigenfunctions are the plane waves 
$(qf;r) = exp(iqfr)/27r in the relative motion with momentum q. The total kinetic energy of 
the two glueballs is 

Esiq) = 6 + esiq) , es{q) = {Q/35) + {2/21) f + 0{k) . (2.54) 

In the heavy particle channel there is no interaction between the bound states. Concerning 
the channel coupling potential, eq. (2.46) carries over as it is since Vcc^ is trivially extended 
to an operator on Tisc- It follows that the role of the free hamiltonian in the space of the two 
channels Tisc ® Ti-hp is played by the operator 

Ho f ^ M . (2.55) 

\^ M y ^ ^ 

Introducing (2.55) the full hamiltonian is split according to 

H8 = Ho + V, (2.56) 

where V represents the interaction among the light particles and the coupling between the 
channels 

(2.57) 

Vcc^ 
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Y = K 



2 



The M0ller operators corresponding to this situation are given by 







in/ out 



Hs, Ho) 



-hm e 



iHst -iUnt 



oi p 



f2.58) 



The reader is reminded that the wave operators are defined by first projecting on the absolutely 
continuous subspace of the free hamiltonian (cf. ref. [17], for example). For the present case this 
is achieved by Pgc, the projection operator on Tisc- Since the continuous subspace is spanned 
by the plane waves the (improper) in- and out-going scattering states are obtained as 







in/ out 




(2.59) 



and it is straightforward to deduce the Lippmann-Schwinger equation. We introduce the Green 
operators 



lim 



1 



G 



1 



hp(gj 



(2.60) 



p\o es{q) - T + ip ' es{q)-M 
that describe the propagation with energy Eg,{q) in the two channels. Then the desired integral 
equation for the in-going scattering solutions of the two-channel system is found to be 





/ 







/ 








(2.611 





Ghp(g) 

If one substitute the second component of this equation in the first one can eliminate the heavy 
particle wave function. Thereby one obtains a modified integral equation for the scattering 
states of the two-particle channel 

$-(g) = $(g) + G,c(g) W(g) $-(g) , (2.62) 

with an energy-dependent potential 



W(g) = V,c + VccGhp(g)V, 



t 4£f 



V.C + V 



add 



(2.63) 



The above result means that, as far as the the two-glueball channel is concerned, the entire effect 
of the coupling to the heavy particle channel is the replacement of the (original) potential Vgc 
by the effective potential (2.63). The physical interpretation of W(qf) is obvious. The additional 
term Vadd describes a transition from Tisc to ?ihp, propagation there according to Ghp(g) with 
energy Eg,{q), and return to Tisc- As a consequence we must be aware of the occurrence of 
resonances in the collision of the light glueballs when the total energy exceeds the mass mi of 
the lowest lying heavy glueball. Furthermore, since the channel coupling is weak we expect a 
possible resonance to be sharp. 

For the scattering matrix one deduces the familiar result 

($,(g'), Uq)) - 2vri/c^ 8 {E,{q') - E,{q)) ($,(g'), W(g)$i"(g) 



S(g',g) 



(2.64) 



where the index s indicates that we did not forget to project on the subspace of symmetric 
wave functions. 
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2.5 Von Neumann series 

We are heading for an approximate solution of the Lippmann-Schwinger equation (2.62) in 
powers of the hopping parameter n. The general idea is the same as discussed in section 3.2 of 
[1]. First we neglect the 0{n) contributions to (2.62) and solve the resulting integral equation 
exactly. This leads to a modified Lippmann-Schwinger equation which is solved by means of a 
von Neumann series in the higher order contributions. 

The leading order channel coupling (2.46) gives rise to a leading order additional potential 

Vidd = VLGhp(g)Vj;. (2.65) 

By means of (2.65) the effective potential (2.63) is decomposed as follows 

W(g) = Wi(g) + 0{k) = VL + Vi,, + 0{k) . (2.66) 

We substitute (2.66) in the Lippmann-Schwinger equation and obtain 

$-(g) = $(g) + GJq) W\q) $-(g) + GJq) [W{q)-W\q)]^ . (2.67) 

Clearly, the third term is of the order 0{k,)^ so we have achieved the desired split of the integral 
equation for the scattering states in SU(2) gauge theory. 

The leading order approximation for the wave function describing the collision of two glue- 
balls is given by the solution of the integral equation 

<^T{q) = Hq) + GsM W^(g) $r(g) . (2.68) 

With this solution at hand the von Neumann series for the scattering states can be derived in 
complete analogy to the discussion in [f]. 



2.6 Leading order scattering solution 

We consider the leading order Lippmann-Schwinger equation in the coordinate space represen- 
tation 

<5Tiq; r) = $(g; r) + ^ GJq; r, r') W\q) <5Tiq; r') . (2.69) 
The Green function has the integral representation 

GJq; r,r')= -P- -r- — -- . (2.70) 

Jb (27r)2 £§(?) - esip) + ip 

It fulfils the identity 

[-T + es{q)] GJq; r,r') = S{r - r') , (2.71) 

and is invariant with respect to planar rotations 0(2, K) in both coordinate arguments. As a 
consequence, the integral equation for the symmetrized scattering solutions is obtained if we 
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replace the wave functions in eq. (2.69) by their symmetric counterparts $™;^(qf; r) and $s(9; ^'), 
respectively. 

We start the evaluation of the RHS of eq. (2.69) by considering the contributions that 
are due to the potential V^^ (cf. eq. (2.53)). The operator Dnc projects on those functions in 
the physical Hilbert space which are zero outside the region |r| < 2. They are very rapidly 
decaying and do not belong to the absolutely continuous spectrum. Hence, the symmetric 
scattering solutions must be orthogonal to these functions, i.e. they vanish at r = k^k ± 1. 
From these considerations we conclude that the operators TDnc, DncTDnc and D(A) do not 
contribute to the Lippmann-Schwinger equation (2.69).^ One finds 

^G',,(g;r,r')VL$i:\(g;r') = - ^ GJq;ri)aq;v) • (2.72) 
r' |77l<2 

The quantities ^[q; rj)^ \ri\ < 2 result from the action of the operators Dt^T on $™;^(qf; r). They 
are sums of values of the symmetric scattering solution at definite coordinate points. Since 
C{Qi v) = C{Qi ~'n) there are seven of these unknown functions in eq. (2.72). 

The additional potential V?^^^ is calculated explicitly by inserting the projection operator 
(2.42) on the heavy particle channel in eq. (2.65) 

VLd = VL G^M PhpVj; = E ^-'^y^-' . (2.73) 

t^i £8{q) - P'k 

The single contributions to this sum are of the type 



|77.)°°(77.| 

9.,k /\ , A; = 1,2,3, 2<|77,|,|77^.|<3, g,^,eM. (2.74) 



For sake of completeness I quote the result explicitly 

VLd = E [(5+|2-U/)° + 5-|2-^-/)°) (°(2-U/|+°(2-/ + ^| 

ki^i=i 

-^(g+\2-k- If + 5_ |2 • A; + if) (°(2 • A; - /| + °(2 • / - ] , (2.75) 

where the coefficients (/_|_ and </_ are given by 

0.000327 0.00572 , 0.00605 
9± = —-^ + —-^ ± • 2.76 

The rather small numbers in (2.76) are due to the weakness of the coupling between the light and 
the heavy particles (cf. eq. (2.47)). It follows that the additional potential will only contribute 
significantly if the energy of the colliding glueballs is near the mass of one of the heavy particles 
|1), |2) or |3). 



^As a check I also did the calculation without using this a priori information about the behaviour of the 
scattering solutions near the origin. Like in the case of the Ising model the result corroborates the prediction, 
I.e. <D-(g,fc) = <D- (g,fc±0 = 0. 
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In the coordinate space representation of the Hilbert space Tigc, (2.74) act as local multi- 
plication operators Dfj, 2 < |77| < 3, possibly multiplied with a difference operator. Like V^^ 
(cf. eq. (2.53)), the effective potential is short range. For notational convenience we introduce 

^^■.(g) = E , Jf , • (2-77) 
k=i ^8{q) - P'k 

and obtain 

Y: GJq; r, r') Vi^a Kiiq'^ r') = Y^g^GsM, r, Vi)Ki(q-^ V,) • (2.78) 

T' i,j 

The above results eq. (2.72) and eq. (2.78) show that the coordinate dependence of the 
leading order scattering states is determined by the free solution and the Green function. As 
unknown quantities there are the sums (^(g;?7), {fjl < 2, and the values ^^^^[q^rj-) of the 
scattering solution in a disc 2 < \ri-\ < 3. Altogether these are II complex valued functions 
which have to be determined in order to solve the leading order Lippmann-Schwinger equation. 
We adopt the notation 

[770, 77i, • • • , ^7io] = [0, i, 2, 2 • i, 2 • 2, i + 2, i - 2, i + 2 • 2, i - 2 • 2, 2 • i + 2, 2 • i - 2] . 

If one substitutes r = rj-^ z = 7, . . . 10, in eq. (2.69) and acts with Dt^.T, z = 0, ... 6, on both 
sides of the Lippmann-Schwinger equation using the identity (2.71), one is left with a set of II 
equations 

6 10 

EM.C(9;^.) + EM.^m(9;^.) = e,{q)^^{q-ri,) for z = 0, . . . , 7 , (2.79) 

1=0 j=7 
6 10 

EM.C(9;^.) + EM.^m(9;^.) = UT^Vi) for z = 7,..., 10. (2.80) 

1=0 j=7 

The elements of the II X II matrix Ai consist of sums of Green functions evaluated at definite 
coordinate arguments, multiplied with the energy-dependent terms gij{q)- Due to its complexity 
an analytical inversion of this object is a formidable task which even the algebraic manipulation 
program MAPLE [18] is not able to do. However, all quantities appearing in Ai are well-behaved 
functions of the relative momentum of the colliding particles. Hence, for any value of q we are 
dealing with a complex II X II matrix that MAPLE is able to invert. This provides us with 
an explicit expression for the unknowns and hence for the leading order scattering solution. 
Introducing 

es{q)M-/ forj=0,...,6 
= I (2.81) 
M-^ forj = 7,...,I0 , 



the result i 



IS 



10 



aq;Vi) = EK,nq;V,) forz = 0,...,6, (2.82) 



10 



^Zi^^Vi) = EK.nq^V,) for z = 7,..., 10. (2.83) 

1=0 
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2.7 Leading order scattering amplitude 

In order to obtain the leading order transition matrix one has to calculate the scalar product 

T,(g',g) = ($,(g'),Wi(g)$i-(g)) =^$^(g';r')W^(g)$i^(g;r') . (2.84) 

r' 

The action of the local effective potential on the leading order scattering solution has just been 
worked out. In fact we only have to replace the Green function GsdQ] f , f') by the symmetric 
free solution $s(qf;r') in eqs. (2.72) and (2.78). One finds 



Ti(g',g) 



-$.(g';0)e(g;0)-2^$^(g';77,)e(g;^. 



8 = 1 



10 



(2.85) 



hi =7 



In view of the forthcoming calculations it is worthwhile to substitute the results (2.82) and 
(2.83) for the leading order scattering solution explicitly. We introduce yet another 11 X 11 
matrix 

-2^f^J for z = 1, 

10 

2 gtkiq) J^kj for z = 7, . . . , 10 

k=7 

In terms of this quantity the result becomes most concise 



.6 



(2.86) 



10 



"n"i(g',g) = ^f^J^s{q'■,r)i)^s{q■,r1j) ■ 

jj=0 



(2.87) 



Since both the Green function (2.70) at definite coordinate points and the terms (2.77) only 
vary with the kinetic energy Eg,{q) of the colliding glueballs the momentum dependence of the 
on-shell leading order transition matrix is completely determined by the free solution $s(9)- 
This remarkable observation considerably simplifies the evaluation of (2.87) and hence the 
diagonalization of the scattering matrix (2.64). 
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3 Phase shifts in lattice gauge theory 

The physical quantities that characterize a scattering process are the phase shifts. Making use 
of the results obtained in the previous sections we can now discuss how the phase shifts for 
the elastic collision of the lightest SU(2) glueballs in two spatial dimensions are determined. 
Naturally, many ideas carry over from the detailed discussion of the (3 + 1) dimensional Ising 
model in [1]. Therefore we can be rather concise from time to time. The exact eigenvalues of 
the scattering matrix are written as exp(2i(5). In order to label the phase shifts properly we first 
consider the low energy regime. At small momentum the rotational symmetry gets restored on 
the lattice and the phase will behave like in the continuum. They can be identified by their 
threshold behaviour which is ruled by the angular momentum quantum number / (cf. appendix 
C.2). For each / one has to distinguish further between the symmetry sectors F of the discrete 
symmetry group in two dimensions 0(2, (cf. appendix C.l). So the phase shifts carry the 
quantum numbers <^/(F). 

3.1 Threshold behaviour 

To begin with we analyse the scattering quantities in the limit where the relative momentum 
of the colliding particles is very small. The kinetic energy was calculated to second order in 
the hopping parameter n 

^2 2 

^8(g) = 2mr + — = 2mr + — + 0(9^) , (3.1) 
mk mk 

where q= \q\. We introduce momentum "circular" coordinates 

= qcos{<p) , q2 = qsin{<p) , 

and define energy eigenstates 

The normalization is chosen such that 

{Ei <f>' I E,, <f>) = {E',, <f>' in I E,, <f> in) = 8{E', - E,) 8^{<f>' - <f>) , 
where we have abbreviated the symmetrized (5-function 

8M' -<P) = (1/2) m -<P) + 6{<P' -<P-7r)] . 
With respect to these bases of states the scattering matrix assumes the following form 

S{E',, <f>'- E,, <f>) = S{Es; <t>\ <t>) 8{E', - E,) , (3.3) 
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with 

S(^8; <t>\ <t>) = W -<t>)- vri^'mk -J{Es; <t>' , <t>) + 0{q) . (3.4) 

For an analysis ol the small momentum behaviour of the leading order transition matrix we 
cannot use the explicit formula (2.87) since the matrix J\f could not be determined analytically. 
Instead we go back to eqs. (2.79) and (2.80). The expressions (2.77) as well as the free solution 
can be expanded in powers of q straightforwardly. As far as the Green function is concerned it 
is shown in appendix B that 



^ ^ 21 



- i + Hv) 

Ztt 



+ 0{q), 6(77) GIR. (3.5) 



Hence, in order to obtain the threshold behaviour of the scattering solution we can neglect 
every term that vanishes like 0{q). In particular we replace all cosines by 1 on the RHS of 
(2.79) and (2.80). Thereby the system of equations considerably simplifies and is easily solved. 
For small relative momentum the leading order scattering solution has the form 



1 

ln{q) + TT c{r ) 



^) = 7-7-^—777 + 0(9) , c(r)GC. (3.6) 



Substituting in (2.85) we find for the leading order transition matrix 

T'(^«^^'*> = ^7Ml3fc^i + OW. (3.7) 

The above result does not come as a surprise since it is the typical threshold behaviour for finite 
range potential scattering in two space dimensions. In order to corroborate this statement I 
discuss a simple quantum mechanical problem in appendix C.2. At low energies the amplitude 
is isotropic in the center of mass system and the scattering is pure s-wave 



27r ln(1.527?) ^ , , ^ , ^ 



For the total cross section one finds (cf. appendix C.2) 

^ = — 2., J,, , +Oiq). (3.9) 

q 41n [1.527 q) + tt^ 

3.2 Energy dependence 

As in the case of the Ising model we take advantage of the fact that the kinetic energy (up 
to the order calculated) is a simple polynomial in the absolute value q = \q\. Concerning the 
leading order transition matrix (2.87) it follows that the most complicated contribution, the 
11x11 matrix A/", only varies with q. Therefore it is sensible to introduce momentum circular 
variables 

q^ = 2 sin((j';^/2) = qcos[(f>) , q2 = 2sin((j'2/2) = qsin[(f>) . (3.10) 
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The associated free and scattering solutions are denoted by (f>) and respectively. 

Their normalization is calculated straightforwardly 



where we have introduced the normalization function 



^^^^ ^ I (1/4)^16- + ?*sin2(^)cos2(^) if - 2 < ?,(?, ^) < 2 

[ otherwise . 

The scattering of the two glueballs is invariant with respect to lattice rotations 0(2, K) 
(this discrete group of planar transformations is discussed in some detail in appendix C.l). 
Consequently, the (symmetric) eigenstates of the scattering operator are grouped together in 
multiplets that transform according to the one-dimensional irreducible representations F = 
Ai, A2, B2 of 0(2, 7L). This decomposition is made manifest by introducing a new basis for 
the free and the scattering states 

$.(?,^) = EE^rW|i^8;r,/) , $in(^,^) = EE^rW|i^8;r,/in) . (3.13) 

r / r / 

The angular functions are defined in appendix C.l. They are a basis for the space of 
symmetric, infinitely differentiable functions on a unit circle. The index / represents the (con- 
tinuum) angular momentum quantum number. For each /, Xf spans a one-dimensional irre- 
ducible subspace for the representation F of 0(2, K). According to the decomposition (C.2), 
the colliding particles carry different total angular momentum in the different symmetry sec- 
tors. Clearly speaking one has / = 0, 4, 8, 12, . . . , oo for F = Ai, / = 4, 8, 12, . . . , oo for T = A2 
and / = 2, 6, 10, . . . , 00 for F = i?2. 

The basis transformations (3.13) reduce the scattering matrix to block diagonal form 



SiEi ct>'- E,, ct>) = 6{E', -E,)Y.Y1 ^ri^') SiE,, F),, X^^) , (3.14) 

r 

where 

S(^8, F),,, = Z{Es, T)i,i - 2tt\k^ T{Es, F),,, . (3.15) 
For the normalization matrix ^(E'sjT) and the transition matrix T(£'8,r) one gets 

Z{E,,T)i,i = 2m^' Jd^Z{E,J)X^{^)Xl{^)\ (3.16) 

T{E„T),i = j A^' j A^xU^')T{E,-A',^)X[{^f . (3.17) 

In order to determine the phase shifts it remains to diagonalize the modified scattering matrix 

S'(^8, r) = Z{Es, F)-i S(^8, F) = ]l - 2TTi^ Z{Es, T)-' T{Es, F) . (3.18) 
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Figure 3.1: The scattering phases in the Ai sector as a function of the momentum q^. 

As long as q < 2 the integration over the angular variable in (3.16) is unrestricted (cf. the 
definition of the normalization function (3.12)), and one can use the orthonormality of the 
angular functions A'f (^) to invert the normalization matrix 

Z{E,,T)j;l = (l/2)mk/ d$Z{E,J)-'Xi:.{$)XU$)* . (3.19) 

The inverse normalization function has a series representation in powers of the momentum 
square q^ . The coefficients are trigonometric functions of the angular variable (/) and can be 
expressed through exponentials. Obviously, the same is true for {(/)). It follows that each 
matrix element (3.19) is a sum of one-dimensional integrals over a product of three exponential 
functions. Those integrals are zero unless the exponents add up to zero. 

Concerning the leading order transition matrix Ti(£'8;^', ^) the angular dependence is 
through the free solutions only (cf. eq. (2.87)). To be explicit, if (?7i,?72) denote the two com- 
ponents of the lattice vector 77, one has 

^s{q; v) = (2vr)"^ cos(?7i5i + 772^2) = (27r)"^[cos(?7i5i) cos{f]2q2) - sin{f]iqi) sin{f]2q2)] ■ 

Now both cos[r]kqk) and sin[r]kqk) can be decomposed into a sum whose contributions are of 
the type 

cos^ (qk) sin" (qk) , l^,v<\f]k\, G INq . 
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Figure 3.2: Partial d-wave phase shifts occur in the Bi and B2 sector. 



From the coordinate transformation (3.10) it follows that 



cos(?fc) = 1 - ?fc/2 , sin(5j = {qj2)^4:- ql . 

Hence Ti(£'8; (p) has a series representation in powers of q as well, and the coefficients are 
trigonometric functions of the primed and unprimed angular degrees of freedom. These terms 
can be expressed as sums of representation functions (^') and Xf (^) such that the integrals 
(3.17) can be evaluated algebraically. If we tolerate some small error we can truncate the series 
for (3.17) and (3.19) at some order 0[q^) (for q = 1.8 we need ~ 50 to guarantee an error 
of less than 1%). As a consequence these matrices become finite dimensional. 

Now the scattering matrix (3.18) can be diagonalized straightforwardly in each symmetry 
sector r. The associated (non- vanishing) phase shifts are plotted in figs. 3.1 and 3.2 as a 
function of the on-shell momentum 



9e 



nT-hiEg — 2mj.) = q + 0{k^ 



(3.20) 



(in [1] the corresponding quantity was named q-,-). These results represent the leading order 
approximation of the scattering phases characterizing the elastic collision between the lightest 
SU(2) glueballs in two space dimensions. We find effects in the partial s-wave, d-wave and even 
g-wave, on the proviso that there may be significant lattice artefacts in the region q^ > 1. The 
observed behaviour indicates a strong attractive interaction between the two glueballs. 
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3.3 The effect of the heavy particle channel 

In the Ai sector, around ~ 1.7, the / = scattering phase shows a behaviour that is typical 
in presence ol a sharp resonance. Indeed, at = 1.707 the two-glueball energy Eg,{q) is equal 
to the mass mi ol the lowest lying excitation in the heavy particle channel. So we are in a 
region where the additional potential (2.75) is relevant. Note that the latter only contributes 
to the Ai sector since the heavy particles (2.38 - 2.40) have the appropriate quantum number. 

Instead of determining the resonance parameters from the data fig. (3.1) one may try to 
investigate the behaviour of the system in the vicinity of the resonance by a suitable approxi- 
mation. If Eg,{q) — nil the dominant contribution to the Lippmann-Schwinger equation comes 
from the additional potential (2.75). Moreover, we can approximate 



9+ - 9- 



0.000327 def 



and obtain 



^ add — 



\2 2 2 

- — E E (12 

/^i k^i=i k'^i'=i 



k + lf + \2-k-lf 
°(2-F + /'|+°(2-F-/'| 



(3.21) 



It is advantageous to consider this quantity as the effective potential of a simplified two-channel 
resonance model (cf. ref. [19]). To this aim we introduce the hamiltonian 




(3.22) 



that acts on the Hilbert space built by the two-glueball configurations Tisc and the heavy state 
|1). The coupling between the two channels is represented by the potential 



Vres = A (|2-A; + /)° + |2-A;-/)°) (1| . 



(3.23) 



Within this system the collision of the light glueballs is determined by the Lippmann-Schwinger 
equation 



/ 







+ 



/ 



V 











I 



V,, 
V ^ 



/ 



(3.24) 



As discussed above (cf. eq. (2.62)) one could easily derive a modified equation for the two- 
particle channel with an effective potential 



V V ^ 



£8(9) 
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that is nothing but the approximate potential (3.21) we started from. 

However, we will step along the complementary path here. We will not eliminate the heavy 
particle wave function but deduce an explicit representation for it. Then, by means of the 
Lippmann-Schwinger equation (3.24) one finds for the on-shell transition matrix 



T(g', q) = S [esiq') - Ssiq)] ($.(g), V^e. ^'"(g) j . (3.25) 

Acting with Vres^ on the first component of eq. (3.24), substituting the second component and 
inserting the result in the second component again, it follows 

^S\Q) ~ P'l ~ Vres vJsc V res 

Now we substitute this expression for the heavy particle wave function in eq. (3.25). The result 
is the famous Breit-Wigner formula describing the rapid energy dependence of the scattering 
amplitude near a resonance (cf. ref. [20], for example) 

The quantities in the denominator are given by 

10 



10 

V{q) = -2\^ J2 ^^^[Gsciq■,V^,VJ) + GscXq■,V^,-VJ)] 

10 . 

= 47rA2 ^ j^A^p8{es{q)-es{p)) ^siq^-nd^sip^-n. 



whereas the numerator reads 

10 

v{q',q) = 4A2 Yl Mq';Vi)Mq;v,] 

= (4/7r)A^ (/o XI {4>') + /, {<!>')) (/o XI (<P) + f, X\^ (<P)] 
with 

/o = V2(2-f,- + |r-^r) , h = -W + i,^. 

In order to determine the parameters of the resonances it remains to locate the pole in eq. (3.27). 
Obviously, since all quantities are known explicitly, this is easily done, and the wanted momen- 
tum value is = 1.712. So the resonance in the collision of the light glueballs occurs at the 
energy 

E^^, = 6 + K^ eres(1.712) = 6 + 0.451 + 0{k^) (3.28) 

and has the width 

Tres = r(1.712) = 0.0148 + 0{k^) . (3.29) 
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Concluding remarks and outlook 



The present paper demonstrates that a systematic strong couphng analysis of glueball scattering 
is feasible. For the case of SU(2) Yang-Mills theory in (2 + 1) dimensions the phase shifts 
describing the elastic collision of the lightest glueballs could be calculated explicitly to leading 
order. A resonance was detected whose parameters could be determined accurately. 

Further investigations might go in two directions. 
The first is to extend the present work to higher orders in the inverse coupling. This may be 
worthwhile since only the leading order approximation of the scattering quantities is known 
so far. Thereby one would achieve further insight into the nature of glueball scattering and 
could clarify the role of the resonance. I do not expect to meet any severe difficulties, although 
the diagonalization of the scattering matrix will be more complicated since the momentum 
dependence of the two-particle energy will probably be less simple to higher orders in the 
hopping parameter. From a technical point of view it would be desirable to develop an algorithm 
that generates the higher order reduced hamiltonians by computer. 

The second and maybe more interesting possibility is the investigation of lattice models that 
are of direct physical relevance. This would mean to go from two to three space dimensions 
(from the point of view of strong coupling perturbation theory the differences between SU(2) 
and SU(3) are negligible). The number of geometrically distinct Wilson loops of length eight 
increases from 5 to f5, and the different classes contain up to 48 states. In addition, the four- 
link excitations receive a further (spin) degree of freedom since the plaquette may be oriented 
in either of the spatial directions. Clearly, the eight-link sector, which is the staring point for 
the derivation of an effective two-glueball hamiltonian, is much more involved than in two space 
dimensions. However, I would like to emphasize that the resulting difficulties are not due to an 
unsuitable method but refiect the complexity of the theory itself. 

From the conceptual point of view all aspects are already present in the two-dimensional 
model. Therefore the procedure carries over almost verbatim and the problem can be tackled 
following the way described in great detail in the present work. In particular, one can be sure 
that the spectrum decouples to higher orders of the hopping parameter such that the effective 
two-particle system simplifies. Still it must be expected that there are a number of resonances 
in the collision of two glueballs. Such a strong coupling analysis would help to clarify an early 
result of Miinster [25] who found that the glueball interaction strength is extremely large. A 
possible explanation could be the existence of a resonance near threshold. 

In principle there is no obstacle to apply the concepts developed in [I] and this work for a 
qualitative investigation of the low-energy scattering of mesons and baryons. One interesting 
question which seems to be within reach of a strong coupling analysis is whether, in the pres- 
ence of light quarks, the decay of glueballs into pions is really suppressed as predicted by the 
phenomenological Zweig rule. 
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Appendix A 



A.l Elementary properties of the group SU(2) 

The Lie group SU(2) consists of all unitary (2 X 2) matrices with determinant equal to 1. The 
corresponding Lie algebra su(2) can be identified with the space of all complex (2 X 2)-matrices 
X which satisfy 

X = -X^, TrX = 0. (A.l) 

The matrix exponential 

exp(X) = ]l + E TT (^-2) 

defines a map from su(2) onto SU(2). The Lie algebra is a real vector space of dimension 3. 
We choose a basis of hermitean matrices T^, a = 1, 2, 3 such that 

TTTaT, = {l/2)6ab . (A.3) 

They satisfy 

[Ta,Tk]=ieabcTc . (A.4) 

where Cabc is the totally antisymmetric symbol on three indices. To be explicit, one often uses 
the familiar Pauli matrices 

Ta = aj2 . (A.5) 
With respect to such a basis any element X G su(2) can be written as 

X = iX^T, , X„ G IR . (A.6) 

For complex functions f{U) on SU(2) a natural inner product can be introduced through 

{f,g)= [ dUfiUygiU) . (A.7) 



Here dU denotes the invariant group (Haar) measure which satisfies d[gUh~^) = dU for all 
g^h ^ SU(2). The associated square integrable functions represent the Hilbert space iv^(SU(2)). 
For any group element g G SU(2) one may define an operator K{g) that translates the arguments 
U of the wave functions 

mg)f]iU) = fig-'U) . (A.8) 

Thereby a representation of the group SU(2) on the linear space iv^(SU(2)) is introduced. Since 
the scalar product (A.7) is based on the invariant integral, the representation is unitary. Using 
the exponential map we define the operator 

d 

dR{X) = — {R[exp(tX)]}^_^ , for all X G su(2) . (A. 9) 
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This equation introduces a representation of the algebra su(2) on iv^(SU(2)). Since K{g) is 
unitary and = Ta the operators 

E„ = dR(r„), a = 1,2,3 (A.IO) 

are hermitean. They are the infinitesimal generators of the transformation (A. 8). Using (A. 9) 
it follows that 

R{g) = exp{iXaEa) if g = exp{iXaTa) , (A.H) 
and one deduces the commutation relations 

[E„,E5] =ie„5,E, . (A.12) 

The operator 

a = l 

commutes with all generators E^ and is hence identified with the quadratic Casimir opera- 
tor of the gauge group SU(2). After choosing a parameterization of the group manifold the 
infinitesimal generators E^ act as differential operators in the group parameters on the wave 
functions / G iv^(SU(2)). In particular, E^ is equal to the Laplace-Beltrami operator up to 
some constant. 

The eigenfunctions of the Casimir operator (A. 13) transform irreducibly under the gauge 
group. Thus, to diagonalize E^ we have to decompose the Hilbert space iv^(SU(2)) into irre- 
ducible subspaces. It is well-known that the irreducible representations of SU(2) can be labeled 
by an integral or half-integral number j, the angular momentum. Eor a given j the represen- 
tation space has dimension n = 2j -\- 1. So there is a one-to-one correspondence between the 
angular momentum and the dimensionality of the representation. Within the context of strong 
coupling perturbation theory it turns out to be more convenient if the irreducible representa- 
tions are labeled by their dimension. Therefore I will adopt this point of view throughout my 
thesis. 

By means of the Peter- Weyl theorem one determines the representation functions 

Dil\u), a,6=l,2,...,n. (A.14) 

Eor any fixed value of the indices n and (a, b) they are elements of the Hilbert space iv^(SU(2)). 
Eortunately we do not need these complicated functions explicitly. What is important is that 
they span an irreducible subspace of iv^(SU(2)) for a definite value of n. Eurthermore the 
Peter- Weyl theorem asserts that there are no further irreducible subspaces in the Hilbert space 
iv^(SU(2)). It is easy to show that the representation functions D^^^ {U) are eigenfunctions of 
the Casimir operator with eigenvalues (l/4)(n^ ~ !)• 
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A. 2 Clebsch-Gordon series 

The irreducible representations of the gauge group SU(2) are classified by the quantum number 
n. On iv^(SU(2)) the representation functions (A. 14) form a basis of the associated irreducible 
subspaces. We introduce a simplified notation 

U:, = D^{U) . (A.15) 

It is often convenient to interpret (A.15) as complex n X n matrices. Special cases are the 
trivial representation n = 1 where the function (A.15) is just equal to 1, the fundamental 
representation n = 2 where (A.15) are SU(2) matrices 

= Ul, , (A.16) 

and the adjoint representation n = 3 where the representation functions can be identified 
with rotation matrices S0(3). By n we denote the irreducible representation that is complex 
conjugate to n 

U:, = UX = U:J . (A.17) 

It is a specialty of the group SU(2) that n and h are equivalent representations, i. e. all matrices 
and are related by the same similarity transformation 

U'^ = CU'''C-\ (A. 18) 

For the case of the fundamental representation one finds explicitly 

^ab = ^bJ = '^aa' ^a'b' '^bb' 5 (A. 19) 

where Cab denotes the antisymmetric symbol on two indices, with £12 = +1. As a consequence 
the trace of SU(2) matrices is real. A special case is the adjoint representation. Here 3 and 3 
are identical representations. 

From a well-known theorem in group theory it follows that the direct product of two irre- 
ducible representations rii and ^2 of SU(2) forms a unitary representation of dimension rii ■ n2 
of the same group. A basis for the product representation rii (x)n2 on the linear space iv^(SU(2)) 
is provided by the set of functions 

Ualb,Klb2 ^ = l,2,...,ni , 02,62 = 1,2,..., ^2 . (A. 20) 

In general, the product representation is reducible. The decomposition of ni(x)n2 into irreducible 
components is described by the Clebsch-Gordon series 

m (g) n2 = (|ni - n2\ + 1) ® (|ni - ^2 1 + 3) ® . . . (ni + ^2 - 3) ® (ni + ^2 - 1) . (A. 21) 

This is just an alternative expression of the well-known result that the addition of two angular 
momenta with eigenvalues ji and j2 gives all j values spaced integrally between |ji — j^l and 
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(ii +i2)- Corresponding to (A. 21) there is a decomposition of the basis functions (A. 20) of the 
product representation 

= E «i! ^2, «2 I ^, «) (n, b\n,,b,;n^,b^) . (A.22) 

n,a,b 

Here the sum extends over all values n that are contained in rii (/) n2 via (A. 21). The Clebsch- 
Gordon coefficients (ni, ai; ^2, 0(2 | can be considered as a ni •^2 X rii ■n2 matrix whose first 
index is the pair (nin2) and second index the pair (na). By convention this matrix is unitary 

E (n,a I ni,ai;n2,a2)(ni,ai;n2,a2 I n',a') = S^n' Saa' , (A. 23) 

ai ,a2 

J2{ni,a^;n^,a^\n,a){n,a\n^,a[;n^,a'^) = S^y^S^^a'^ . (A. 24) 

Throughout this thesis we only need some of the Clebsch-Gordon coefficients explicitly. They 
can be determined by means of a general theorem which states that the product of an irreducible 
representation and its complex conjugate contains the trivial representation exactly once and 
the adjoint representation at least once in its Clebsch-Gordon series (cf. ref. [21, sect. 16.6], for 
example). For the case of SU(2) the associated coefficients are given by 

(n, a; n, 6 I 1) = —= S^b 

'n 



and ^ ^ 

(2,a;2,6 I 3,j) = -^(cTjOfea , (3, a; 3, 6 | 3, j) = ej„6 , etc. 

Making use of the equivalence relation (A. 19) for the fundamental representations 2 and 2 one 
deduces 

(2,a;2,6 I 1) = (2,a;2,6 | 1) = -^e^^ , 

(2, a; 2, 6 I 3, J ) = -^{aj)bceca , 



(2,a;2, 6 I 3, j) = —ebc{crj)ca- 



1 

7!' 



A. 3 Strong coupling matrix elements 



The matrix elements in strong coupling perturbation are such that a product of several link 
operators corresponding to irreducible representation functions is matched between the static 
vacuum 

I K\b, Klb2 ■ ■ ■ Klb, I ^) ■ (A.25) 

These objects are zero unless the decomposition of the product contains the trivial represen- 
tation. The actual value of (A.25) is calculated straightforwardly using the Clebsch-Gordon 
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series (A. 21) and (A. 22). In the following I give a list of those matrix elements that contribute 
to the first and second order processes relevant throughout this thesis. 
The simplest case is the product of two fundamental representations 

2® 2 = 1 ®3 , 

and the associated non- vanishing matrix element is 

(0 I \]ah\h'h' I 0) = (1/2) Saa' Sw • (A.26) 

In case that the complex conjugate representation 2 occurs instead, the corresponding matrix 
element can be easily deduced from (A.26) by means of the similarity transformation (A. 19), 
for example 

(0 I l^ah^a'h'^ I 0) = (1/2) 8^' • (A.27) 

Three times the fundamental representation does not contain the trivial one, so there are no 
contributions from such matrix elements. In second order processes we encounter the product 
of four representation operators. The Clebsch-Gordon series tells us that the intermediate state 
may be in the trivial or adjoint representation 

2®2®2®2 = (1®3)®(1®3) = 1®1®3®3®3®5. 

Let P(l) and P(3) denote projection operators on the corresponding irreducible subspaces. The 
associated matrix elements are 

(0 I Ua6 Ucd P(l) Ua'6' Uc'd' | 0) = {I I ^) £ ac ^bd £ a' c' ^b' d' , 

(0 I Ua6 Ucd P(3) Ua'6' Uc'd' | 0) = (1/12) {2 e ac' £ a' c — £ ac £ a' c') 

X {2ebd' £b'd — £bd£b'd') ■ 

It is only for the investigation of the seven-link one-particle subspace that more complicated 
contributions occur. In particular one needs 

I U., U^, P(2) U,,,, U^,,, I 0) = (1/6) 6^^, 6,,, 6^, 6,,, , 

I ^ab Ucd P(4) U^,;,, U^,^, I 0) = (1/12) e^f^ e,,,^ e^^,^ ((^^)a'a • 
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Appendix B 



The Green function (2.70) determines the coordinate dependence of the scattering solution 
describing the collision of two SU(2) glueballs in two dimensions. In order to calculate the 
corresponding transition matrix it has to be evaluated at various coordinate points. Substi- 
tuting the second order approximation of the kinetic energy (3.1) one obtains the integral 
representation 

G4,;r) = (21/2)/ ^.1^41^. (B.l) 

Mind that this is the two-dimensional counterpart of [1, eq. C.3], so we can adopt many argu- 
ments from appendix C of [1]. But cautiously, the dimension is different. 
Clearly, there is a representation in terms of Bessel functions 



GUq.r) = -{21 1 A)\ AtFUt), (B.2) 

with 

F,c(t) = exp [(i/2) ((t - 4) t] e-"' H ^''"JrAt) ■ (B.3) 

k=i 

Be careful not to interchange the t-integration and the limit — )• here because the asymptotic 
behaviour of a product of two Bessel functions is . In detail, the asymptotic expansion reads 

2 oo 

^-(t) = E E K^i^) exp[-(^, + p) t] t-^-' , z, = (i/2) (8 - 4/i - q') . (B.4) 

^=0 i/=0 

The series that is truncated after u = N is denoted by F^{t). As described in detail for the 
three-dimensional case we estimate the rest (cf. [1, eq. (C.ll)]) and obtain an approximation 
for the Green function with arbitrarily small error e 

G',,(g;r) = -(21/4)i \ T dt Fjt) + T dt F,^(t) 1 + (9(e) . (B.5) 

To be explicit I demand e = 10~^^ and choose N = 20. Then the limit of integration turns out 
to be tjv = 17 for r = 0, for example. 

As far as the evaluation of the finite integral is concerned we can forget about the limit 
— )• in (B.3), and the reader is referred to appendix C of [1] once more. We turn to 
the infinite integral in eq. (B.5). Unless p = the integration over the summands of the 
asymptotic expansion (B.4) is well-defined even without the convergence factor e"''*, and we 
obtain a formula similar to [1, eq. C.13]. So it is only the contributions oc from the 
asymptotic expansion of the integrand that make it necessary to discuss the calculation of the 
Green function in two space dimensions separately. By means of ref. [22, sections 3.381, 8.214] 
one gets 

9 

^ /'CO 

lim V b^r) / dt exp [-{z^ + p) t] 

,,-n JtN 



11=0 

2 2 

li^^ E bi,o{r)Fi [-(^^ + /?) tjv] = - E bi,o{r] 
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Hz,t^) + c + Y^ 



.=1 k-kl 



where C = 0.577... is Euler's constant. For computational purposes the infinite sum is truncated 
such that the error is less than e. 

The above result shows that the Green function (B.l) is logarithmically divergent if = 0, 
that is for (f^ = 0,4 and 8. In particular, one finds the threshold behaviour 



GsA,q;r) = (21/8) 



27r 



i + Hr) 



+ 0{q) , b{r) G IR 



(B.6) 



Contrary to the Ising model, where the Green function was needed only at r = and 
r = 1 for our purposes, the calculation of the transition matrix (2.87) affords the evaluation of 
Gsc{q] f) at a considerably larger number of spatial points (cf. eqs. (2.72) and (2.78)). Although 
the method described above allows to determine the Green function for any value of the relative 
momentum q and the lattice vector r it is much more elegant and economic to complete this 
straightforward calculation by some algebraic considerations [24]. 

Consider the integral representation (B.l). First of all it is easy to deduce the identity 
(cf. eq. (2.71)) 

(2/21) i-A + q') G.,{q;r) = 6{r) . (B.7) 



fB.8) 



By means of a partial integration it is equally straightforward to show that 



Gsciq; r + k) - Gsciq; r - k) = Vk H{q; r) 



A; = 1,2 , 



where H[q;r) does not depend on the spatial component Vk- Joining (B.7) and (B.8) together 
and assuming that r ^ we can determine the unspecified function 



H{q;r) 



1 



ri + r2 



4-f) G,,{q;r)-2j2G,,{q;r-k) 



k=i 



(B.9) 



provided that r\ ^ ^ 0. After substituting this result in (B.8) again one is left with a 
recurrence relation for the Green function 



G^^^q; r + k) = Gsc{q; r - k) + 



ri + r2 



4-^2) GUq;r)-2Y,GUq;r-k) 



k=i 



(B.IO) 



Due to the planar symmetry 0(2, 7L) we can restrict attention to the sector ri > r2 > 0. Within 
this quadrant eq. (B.IO) allows one to express GscXt^'"') ^ linear combination of the values 
at the corners of the unit square 



ac(g;0), G^q-l), ac(g;i + 2) . 

Furthermore, assuming that r = one deduces from (B.7) 

ac(g; i) = (1/4) U - f) G^^iq- 0) + (21/8) 



:B.121 
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Appendix C 



C.l The rotation group in two dimensions 

The spatial symmetry group of the continuum (2 + l)-dimensional gauge theory is the group 
0(2). In addition to the proper rotations about the z axis it contains reflections in a plane 
containing the z axis. Therefore this group is non-abelian. The irreducible representations can 
be labeled by positive integers / = 0, 1, . . . , oo which are interpreted as the quantum numbers 
representing the (spatial) angular momentum. Except for the trivial representation / = 0, 
which is one-dimensional, all others are of dimension two. The space of infinitely differentiable 
functions defined on the unit circle is an infinite-dimensional homogeneous space for the action 
of the group 0(2). The 'circular harmonics' 

f 1/^2^ for / = , m = 

I exp(im/^) I \/2tt for / =I,2,...,oo, m = ±1 

span the irreducible subspaces which we denote by T;. 

On a two-dimensional spatial lattice the symmetry group is reduced to the planar rotations 
0(2, K). This point group is the symmetry group of the square with undirected sides and 
is often called the dihedral group 1)4 (cf. ref. [21], for instance). It is a discrete group of 
order 8 with 5 irreducible representations Ai, A2, B2 and E of dimensionality 1,1,1,1 and 
2. Through the obvious restriction the irreducible representations / of the continuous group 
0(2) are also representations of the point group 0(2, K). The corresponding representation 
functions are the same modulo 4, i.e. 



DiZ^\R) = D^nl'miR) for all n G IN, i? G 0(2, 7L). 

The decomposition of each of these representations in terms of the irreducible representations 
of 0(2, K) is 

= Ai 

1 = E 

2 = 5i®52 

?, = E 

4 = Ai ® A2 _ 

For lattice calculations it is useful to introduce new basis functions for the subspaces T; which 
make the above reduction manifest. In the context of two-glueball scattering we can restrict 
ourselves to even parity solutions, i.e. even / values, because the wave functions are symmetric. 
As a consequence the E'-sector, which only contributes to odd subspaces T;, is irrelevant. We 
define (cf. ref. 



(C.2) 
> modulo 4 . 
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= cos(/<^)/0F 

X' ((^) = -isin{l<P)/y^ J ' ' ' (C.3) 



for / = 4,8,12,... 



X's^{<P) = cos(/<^)/0F 
XU<I>) = -i sin(/<^)/V^ 



for / = 2,6,10,... 



These functions are a basis for the space of symmetric, infinitely differentiable functions on a 
unit circle 

r / 

where the sums run over the representations F = Ai , A2 , i?i , i?2 and the associated angular 
momentum values. Furthermore, they are orthonormal 



C.2 A simple quantum mechanical example 

One or the other of the readers may be not so familiar with quantum mechanical potential 
scattering in two space dimensions. Those are invited to go through a simple example which 
reveals some of the kinematical peculiarities. 

Consider the relative motion of two identical, non-relativistic bosons of mass m interacting 
via a rotationally invariant potential 

( -Vo < for r <R 

V(r)= (C.4) 
[ for r > i? 

Analogous the well-known three-dimensional case (cf. ref. [23], for instance) the partial wave 
expansion for the wave functions reads 



Hr) = j:i^iJr)Yr{<P) . (C.5) 

Im 

The radial components can be written as 

4^im{r) = Q„ ui{q; r) , (C.6) 

where is a suitable normalization factor and ui[q; r) solves the radial Schrodinger equation. 
The latter has two linearly independent solution given by the Bessel functions J/(r) and the 
Neumann functions A^/(r), for example. 

The physical wave functions are uniquely determined by means of the regularity condition 

lim r~' ui(q;r) = 1 . (C7) 
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In case there is no interaction [Vq = 0) this condition selects the Bessel functions as the regular 
solutions, that is the wave functions describing the independent motion of the two particles. 
For large r, outside the interaction range, the scattering solutions must behave like the free 
ones, hence 

/ X COS [qr - (//2)7r - 7r/4 + 6i{q)] 
M<i;r) Ci{q) -= . (C.8) 



Here Ci{q) is some coefficient, and the phase shifts 6i[q) reflect the effect of the potential 
compared to the free solution. 

Now our demonstrative problem can be solved easily. Outside the range of the potential 
the regular solution is a linear combination of Bessel and Neumann functions 

ui{q;r) = ai{q) Ji{qr) + I3i{q) Ni{qr) for r > i? . (C.9) 

Substituting the asymptotic expansion for Ji{qr) and Ni{qr) (see ref. [22, section 8.451]) and 
comparing with (C.8) one obtains a relation between the coefficients ai{q) and fii{q) and the 
phase shifts 6i{q) of the system 

lc.■c( ^^ M<l) - i A(g) /p-in^ 

exp[2i6i{q)] = ... or ^ ■ 
ai{q) + ipi{q) 

Inside the interaction region the scattering solution is 

ui{q;r) = ^i{q) Ji{\r) for r < i? (CH) 

with \^ = q^ -\- mVo. Demanding continuity for the wave function and its first derivative at 
r = the unknown coefficients in (C.9) can be expressed in terms of (C.II). Substituting 
these results in (CIO) the eigenvalues of the scattering matrix are found to be 

exp [2i 8 iq)] = - 'iMmHl'^'W -^J'i{mH^\qR) 
qM\R)H\'\qR)-\Jl{\R)H\'\qR)' 

where we have introduce the Hankel functions 



i7(^)(r) = J,(r) + iiV,(r) , H^\r) = J i{r) - \ N i{r) . (C.I3) 

The prime denotes the derivative with respect to r. 

We analyse the threshold behaviour. For small relative momentum one can replace A by 
Ao = \/mVo and make use of the series representations for the Bessel and Neumann functions 
(cf. ref. [22, sections 8.440, 8.403]) 



TT No{z) = 2,h{z) [Hz 12) + C] + 0{z 



2^ 



^Ni{z) = 2Mz) [Hz 12) + C] - E z^'-' + 0{z') for/>0 
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Substituting in (C.12) the scattering phases become 

27r ln(i?o q) 



wi 



anc 



ith 



Ro = exp 



41n2(i?o9) + 7r2 
Jo(Aoi?) 



+ 0{q) 



C + 



XoRJi{XoR) 

I Ji{XoR) — Xq R J'i{XqR) 
2'(/- l)!r(/ + 1) ' UliXoR) + XoRJIiXoR) 



R 

'2 ' 



TT [Rq 



v2/ 



+ 0{q 



2l+2\ 



for / > . 



As for the well-known three-dimensional case, where the threshold behaviour of the phase shifts 
is Si{q) oc q'^'~^^, s-wave scattering is dominant at low energies. So it is only the / = channel 
that contributes to the scattering amplitude near threshold 



T{q-<l>', 



-m 



27r ln(i?o q) — ivr' 



+ Oiq) . 



(C.14) 



For geometrical reasons the cross section in two spatial dimensions receives a factor q ^ (see 
ref. [20, section 8], for example) 



d 



(7 



d. 



(27r)3 m2 „ 
\T{q;<P\ 

q 4 



Substituting (C.14) the threshold behaviour of the total cross section is found to be 

47r2 1 



(7 



q 4 1n'(i?o9) + 7r^ 



+ 0{q) 
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